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CHAPTER  1 
INTRODUCTION 


Review  of  Phase  I  Studies 

This  report  describes  studies  conducted  for  Phase  n  of  the  project  ‘Restraining  Post- 
Lique&ction  Flow  Deformations”.  To  provide  a  context  for  this  report,  the  background  of  the 
project  will  be  reviewed  briefly  and  some  key  findings  from  the  first  phase  will  be  presented. 

In  Phase  I  of  the  project,  studies  were  focussed  on  the  use  of  piles  to  restrain  flow 
deformations.  This  topic  was  of  considerable  interest  to  the  US  Army  Corps  of  Engineers 
because  pile-nailing  of  the  upstream  slope  of  Sardis  Dam  in  Mistisappi  was  being  conadered  as 
an  option  for  restraining  sliding  of  the  slope  on  a  potentially  liquefiable  thin  layer  in  the  foundation 
(Finn  et  al.,  1991). 

The  key  factors  controlling  the  feasibility  and  cost  of  pile  installations  to  restrain  flow 
deformations  are  pile  length  and  spacing,  stiffliess  and  strength  of  unliquefied  soils  surrounding 
the  piles,  residual  strength  of  liquefied  soils,  the  geometry  of  the  structure,  and  the  intensity  of 
shaking  after  liquefaction  has  occurred.  The  ability  to  analyze  such  a  complex  problem  vdiile 
taking  into  account  nonlinear  behaviour  of  soil,  potentially  large  strain  deformations  in 
unremediated  parts  of  the  structure  and  a  realistic  interaction  between  piles  and  soil  during  both 
static  and  seismic  loading  is  the  essential  requirmnent  for  determining  the  best  location  for  the 
piles,  an  iq)propriate  length  and  tize  and  for  categorizing  the  efiects  of  soil  properties. 

Very  little  is  known  about  the  behaviour  of  pUes  under  these  complex  conditions.  Most  of 
the  evidence  is  from  Japan  where  pile  foundations  have  been  severdy  damaged  in  liquefied  ground 
as  a  result  of  ground  displacements.  However  these  piles  were  detigned  for  vertical  static  loading 
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only  and  both  piles  and  the  coimections  to  the  i^e  ci^s  were  inadequatdy  designed  to  resist 
horizontal  loading.  Piles  in  Oakland  Harbor  were  nmilarly  damaged  during  the  Lonu  Prieta 
Earthquake  of  1989.  Therefore  there  is  a  need  f(Mr  detailed  aiudytical  studies  of  pUe  installations 
to  provide  much  needed  information  on  the  potential  performance  of  piles  undo*  these  demanding 
situations  and  hence  provide  a  framework  of  understanding  for  design  of  cost  effective 
remediation  measures. 

The  computer  programs  TARA-3  (Firm  et  al.,  1986)  and  TARA-3FL  (Finn  and 
Yogendrakumar,  1989)  were  used  to  conduct  the  studies  in  Phase  I.  It  was  the  first  time  such 
arudyses  had  been  performed.  The  findings  of  the  Phase  I  studies  were  used  to  develop  the  dengn 
requirements  for  the  piles  for  remediating  Sardis  Dam.  The  pile  design  resulted  in  very  substantial 
savings  in  remediation  costs  compared  to  alternative  proposals.  The  design  would  not  have  been 
feasible  without  the  Phase  I  studies. 

The  studies  for  Sardis  Dam  could  be  conducted  using  the  plane  strain  2-D  arudyses  in  the 
TARA-3  suite  of  programs  because  remediation  was  necessary  across  a  long  longitudinal  section 
of  the  dam.  Where  such  conditions  do  not  hold,  3-D  analyses  are  necessary. 

Full  3-D  nonlinear  dynamic  aiudysis  of  pile  groups  is  not  a  feasible  proposition  for 
engineering  practice  at  present.  It  makes  impractical  demands  on  computational  speed  and 
capacity  and  renders  the  extension  to  dynamic  effective  stress  analyas  extreme^  difficuh.  Full 
3-D  arudysis  also  irdiibits  the  detailed  parametric  studies  udiich  are  so  important  in  exploring  cost- 
efficiem  remediation  options.  A  simplified  model  for  the  3-D  soil  continuum  under  horizontal 
shaking  by  vertically  propagating  shear  waves  has  been  developed  vtdiich  overcomes  these 
difBculties.  This  simplified  3-D  model  captures  the  significant  motions  and  stresses  in  fijundation 
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soils  with  accq)table  accuracy  during  ground  diaking  and  provides  the  basis  finr  significant 
advances  in  the  sdsmic  analyas  of  piles. 

This  report  is  restricted  to  describing  the  evolution  and  validation  of  a  proposed  8inq)lified 
3-D  method  of  analysis  of  the  dynamic  response  of  piles  that  greatly  reduces  the  conquitational 
demands  for  solving  practical  problems.  It  is  not  a  state-of-the-art  study  of  the  dynamic  reqmnse 
of  pile  foundations.  References  to  the  literature  are  restricted  primarily  to  case  studies  used  for 
validation. 

Structure  of  Report 

A  general  simplified  model  of  the  dynanuc  response  of  the  soil  continuum  is  presented 
first.  Then  particular  models  based  on  different  ^plifying  assumptions  (Matsuo  and  Ohara, 
1960;  Veletsos  and  Younan,  1994;  Finn  and  Wu,  1994;  Finn  et  al.,  1994a,b)  are  derived  from  the 
general  case.  The  performance  of  these  models  is  assessed  using  the  dasric  solution  of  Wood 
(1973)  for  the  dynamic  pressures  against  rigid  walls  by  a  homogeneous  elastic  badcfill.  The 
validation  study  indicates  that  the  proposed  simplified  model  performs  very  wdl  over  a  wide 
range  of  soil  properties.  The  solution  for  rigid  walls  has  important  applications  in  remediation 
studies  in  its  own  right.  A  generally  useful  remediation  measure  is  the  indurion  of  concrete  plugs 
such  as  slurry  walls  to  restrun  sofl  deformations.  Assessment  of  the  stability  of  sudi  walls 
requires  estimates  of  the  seismic  lateral  pressures.  This  type  of  remediation  was  one  of  the 
options  considered  for  Sardis  Dam. 

The  model  is  modified  to  accommodate  piles  in  the  continuum.  The  equations  of  motion 
for  the  soil-pile  system  are  formulated  in  terms  of  firute  demoits.  Elastic  solutions  for  pile 
impedances  are  validated  against  the  full  3-D  solutions  by  Kaynia  and  Kausd  (1982). 
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The  model  is  extended  to  nonlinear  response  by  ensuring  conqMtibility  betwem  soil  strains 
and  the  strain  dq>endent  nxxluli  and  damping  of  soils  continutMisly  during  dynamic  analyw.  This 
is  a  nuyor  modification  of  the  equivalent  linear  method  used  in  computer  prognuns  such  as 
SHAKE  (Schnabel  et  al..  1972). 

The  proposed  method  of  analysis  is  also  validated  using  data  firom  a  forced  vibration  test 
on  a  Frank!  pile  conducted  at  the  Pile  Research  Facility  of  the  University  of  British  Columbia  (Sy 
and  Siu,  1992).  Data  fiom  centrifuge  tests  on  pile  foundations  under  strong  shaking  conducted  at 
the  California  Institute  of  Technology  (Gohl,  1991;  Rnn  and  Gohl,  1987)  allow  validation  of  the 
proposed  method  of  analysis  vdien  nonlinear  soil  effects  are  important.  For  the  first  time,  the 
variation  of  pile  stififiiess  and  damping  with  time  are  traced  during  strong  shaking. 

Finally,  extensions  of  the  method  will  be  discussed.  These  extensio.is  include  analysis  of 
pile  groups,  rocking  effects  and  dynamic  effective  stress  analysis. 
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CHAPTER! 

SIMPLIFIED  EQUATIONS  FOR  DYNAMIC  RESPONSE  OF  FOUNDATION  SOILS 


The  soil  is  assumed  to  be  a  homogeneous,  isotropic,  dastic  solid  with  a  shear  modulus  G 
and  a  Poisson’s  ratio  v.  The  equation  of  motion  for  the  soil  continuum  in  the  horizontal  direction, 
X,  is  written  as 


9cJx 


>Kdiere  o.  is  the  normal  stress  in  the  x  direction  and  is  the  shear  stress  in  the  x-y  plane. 

For  two-dimensional  plane  strain  conditions,  the  stress  components  are  related  to  the 
displacements  by 


2G  du  dv 

=7'^[0~v)— +v— ] 
X  l-2v  '  'Sx  dy 


2G  dv  Ai 

y  i-2v  fly  dK 


flu  dv 


(2.4) 
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Then  the  governing  equation,  Eq.  2.1,  fcH*  undamped  fiee  vibraticMi  of  the  backfill  be  written 


as 


d^u  d^u  d^u 
0G — 5"+G — T  =  P — T" 

ay*  at* 


(2.5) 


^ere  p  is  the  mass  density  of  the  backfill  soQ,  t  is  time  and  6  is  a  fiinction  of  Poisscm’s  ratio  v. 

The  expression  for  8  depends  on  the  iq)proximations  used  in  moddling  the  dynamic 
re^nse  of  the  soil  continuum.  Three  different  approximations  will  be  considered. 

i)  v*0 

It  is  assumed  that  there  is  no  (vertical)  displacement  in  the  y-direction.  Applying  this 
assumption  to  Eq.  2.2  and  Eq.  2.4  gives 


2(l-v) 

i_2v 


(2.6) 


T 


xy 


(2.7) 


Substituting  Eq.  2.6  and  Eq.  2.7  into  Eq.  2.1  and  comparing  Eq.  2.5,  one  finds 

.  2(l-v) 


l-2v 


(2.8) 
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u)  Oy  »  0 

Here  h  is  assumed  that  the  dynamic  normal  stresses,  Oy,  in  the  vertical  direction  are 
negligible.  .^)plying  the  assumption  to  Eq.  2.3,  one  finds  that 


dv  V  fid 
dy  1-vftt 


(2.9) 


Then  fi’om  Eq.  2.2  and  Eq.  2.4  one  obtains. 


2  ^au 

Ov  = G — 

1-v  dK 


(2.10) 


^xy  ^  a^U  ^  V  a^U 

dy  dy^  l-vay2 


(2.11) 


Substituting  Eq.  2. 10  and  Eq.  2. 1 1  into  Eq.  2. 1  and  comparing  with  Eq.  2.S,  one  finds 


0  = 


2-v 

l-2v 


(2.12) 


iii)  Proposed  Method  (Wu,  1993) 

In  this  method  the  condition  ay  -  0  is  combined  with  the  assumption  that  the  layered  half¬ 
space  behaves  as  a  shear  beam.  The  latter  assumption  implies  that 
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X 


xy 


(2.13) 


The  normal  stress  o.  is  found  by  assuming  Oy  *  0  in  Eq.  2.3. 


ax  = 


2  ^au 

1-v^ax 


(2.14) 


Substituting  Eq.  2. 13  and  Eq.  2. 14  into  Eq.  2. 1  and  comparing  with  Eq.  2.5,  one  finds 


(2.15) 


These  three  different  approximations  to  the  response  of  the  layered  half-^ace  are  incorporated  in 
the  general  expression  for  the  equations  of  motion  (Eq.  2.5)  by  the  parameter  6. 

The  capability  of  different  approximations  to  the  response  of  the  soil  continuum  will  now 
be  verified  u«ng  solutions  to  the  classic  problem  of  dynanuc  pressures  against  a  rigid  wall. 
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CHAFTERS 

DYNAMIC  ANALYSIS  OF  RIGID  WAU^SOIL  SYSTEM 


Figure  3.1(a)  shows  the  geometry  of  the  problem  and  its  boundary  conditions.  A  uniform 
elastic  soil  layer  is  confined  by  two  vertical  ri^d  walls  at  its  two  side  boundaries  and  a  rigid  base. 
The  soU  layer  has  a  total  length  of  2L  and  height  of  H.  The  original  wall-soil  problem  can  be 
represented  by  half  the  structure  because  of  the  antisymmetric  conditions.  The  equivalent  problem 
is  shown  in  Fig.  3.1(b),  and  this  is  the  phyacal  model  that  vdll  be  analyzed.  The  ground 
acceleration,  Uq  (t)  is  input  at  the  base  of  the  wall-soil  system. 

Assume  the  displacement  solution  of  Eq.  2.5  has  the  form 

u(x,y,t)  =  SZ  (A  sinanjX+B  cosan,xXC  fflnbay+Dcosbny)  Yn„(t)  (3.1) 

.^plying  the  boundary  conditions 

u  =  0 . aty  =  0 

u  =  0 . atx  =  0 

The  constants  B  and  D  are  determined  to  be  zero,  so 

u(x,y,t)  =  XZ  C,  sinaniX  sinbny  Yoa(t)  (3.3) 


(3.2a) 

(3.2b) 


dK 


=EZ  c,.  „  •  cosEgiX  -  sin  bjiY  •  (t) 


~  =  EZ  Cib„  sina„x  cosb„y.Y^(t) 


Applying  the  other  two  boundary  conditions 


So 

—  =  0 . at  x  =  L 

Sx 


and 


Su 

^=0 . at  y  =  H 

dy 


one  obtains 


and 


an.-cosa„L  =  0 


b„  •cosb„H  =  0 


therefore. 


(2ni-l)x 


2L 


and 


(2n-  l)x 


2H 
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(3.4) 

(3.5) 

(3.6a) 

(3.6b) 

(3.7a) 

(3.7b) 

(3.8a) 

(3.8b) 
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The  mode  shape  functions  are  written  as 

®mn(x.y)  =  CiMna„xsinbay  (3.9) 

and  the  displacement  solution  becomes 

u(x,y,t)  =  ZZ  ®nm(x.y)  •  Ynm(0  (310) 

Substituting  Eq.  3.10  into  Eq.  2.5,  one  obtains 

-G(0aj,  +bj).  Y^(t)  =  pY,^(t)  (3.11a) 

2  ,  u2\  ^inn(0  2  « 

(0®m  ~  V  (3.11b) 

P  Ynm(t) 

The  natural  frequencies  of  the  system  are  found  to  be 

oL=^(0am+l>5)  (312) 


The  frequency  of  the  first  mode  is 


13 


^  n.n"  >> 


(3.13) 


In  the  case  of  undamped  forced  vibn^on  caused  by  a  ground  aocderadon  Uo(t),  the 
governing  equation  becomes 


5^u  5^u  5^u 

P^-(eG^+G-Y)=-puo(t) 
dr  diT  0y^ 


(3.14) 


Substituting  £q.  (3.10)  into  Eq.  (3.14),  multiplying  the  equation  by  the  mode  sluqm  fimctions,  and 
integrating  over  the  domain,  one  obtiuns 


/J  pZ2:<I>«t^(t)  <I>™,dxdy+JJ  EZ  G(0a?+b?)®^Y^.4)^dxdy 

=  -Uo(t)ff  pl>nm(x»y)  dxdy 


(3.15) 


Applying  the  orthogonality  conditions  and  recalling  Eq.  (3.12),  one  obtains 


JJ  P<*i»<*xdy'Ynm(t)+J/  P<l4n<*X<fy  ®mYnm(t) 
=-uo(oJJ  pi>iiin(x.y)*dxdy 


(3.16) 


Yna(t)+(D5m  •Yna(t)  =  -ao(t)  a, 


(3.17) 
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where 


®I1U1  ~ 


Jf  psin(an,x)  sin(b„y)dx  dy 

II  psin^(a„x)-8in^(bay)dx-dy 


®^inn  ~ 
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(2m-lX2n-lK 


Let  ¥■.(!)  =  a«,  •  €*(0, 

fnm(0  +®L  •  fmnO)  =  ““oCO 

For  a  damped  forced  vibration  of  the  wall-soil  system,  a  constant  modal  damping 
introduced  (Seed  and  Idriss,  1967) 


fnin  (0 + 2Xo  nm  •  fim  (0 +®im  •  finn  W  = -fio  (0 

For  a  given  groimd  exdtation  uo(t),  a  close-form  solution  of  the  system  is  found  to  be 

u(x,y,t)  =  SZ  8ina„x-sinbny-aani-fn„(t) 


(3.18) 

(3.19) 

ratio  X  is 

(3.20) 

(3.21) 


vdiere  ^(t)  is  the  time  history  solution  of  Eq.  (3.20)  corres^nding  to  a  particular  modal 
frequency  (On-  It  is  noted  that  Eq.  (3.20)  is  the  standard  damped  vibration  equation  of  a  single 
degree  of  freedom  system. 
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The  dynamic  earth  pressure  acting  on  the  wall  is  the  normal  stress  at  W).  The  dynamic 
pressure  distribution  along  the  wall  is 


dU| 

p(x.y.t)x-o  =  PG^|x-o 

p(x,y,t)x.o  =  PGXZ  *ni®inii  •sin(bny)fmn(0 


(3.22) 


From  Eq.  2.6,  Eq.  2. 10  and  Eq.  2. 14,  the  stress  coefficient  ^  for  the  different  sod  models  is 

Forv  =  0;P»2(l-v)/(l-2v) 

ForOy  =  0;3  =  2/(l-v) 

For  the  proposed  method;  P  =  2/(l-v) 


The  total  dynamic  thrust  acting  on  the  waU  is 


fH 

P(t)  =  i  P(x,y,t)x^dy 


P(t)=PG  lI 

“n 


P(t)=PGSZ 


16fn.(l) 

(**(2ii-1)*)L/h' 


(3.23) 


The  total  dynamic  moment  acting  at  the  base  of  the  wall  is 
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M(t)  =  r  y  p(x,y.t)x^  dy 
0 


M(t)  =  pGZZ 


•a 


bi 


(3.24) 


For  a  harmonic  input  Uo(t)  *  A««  *  e’**,  the  amplitude  of  the  steady  state  reqxmse  ^t)  is  found 
from  £q.  (3.20)  to  be 


f _ ^max 

““  (caJm-fl>*)+2iX(D 


ftm 


‘(D 


(3.25) 


For  any  excitation  uq  (t)  the  time  history  of  the  modal  dynamic  thrust  associated  with  a  particular 
mode  is  obtained  uang  Eq.  (3.23).  The  time-lustory  of  the  dynamic  thrust  for  the  desired  number 
of  modes  is  then  determined  by  summation  of  the  time  histories  of  the  nurdal  dynamic  thrusts. 

For  earthquake  motion  the  peak  modal  thrust  acting  on  the  wall  assodated  with  a 
particular  frequency  coai  can  be  determined  using  the  pseudo-spectral  velodty  Sv***.  The  pseudo- 
spectral  velocity  Sv**  is  derived  from  response  spectral  di^lacement  Sd""  by 

Sr=®n«Sr  (3.26) 


where  S/"  is  also  the  peak  CM(t)  corresponding  to  an  excitation  frequency  <Om»-  From  Eq.  3.23, 
the  peak  modal  thrust  P«  is  determined  by 


17 


=  PG 


(3.27) 


The  peak  dynanuc  thnist  is  estimated  by  oMBbiniqg  the  individiial  peak  modal  tlvusts  by 
some  q)proximate  method.  Either  the  absolute  summation  or  the  root  square  summatkm  of  the 
peak  modal  thrusts  is  commonly  used. 

Static  Ig  Solution:  Validation  of  Model 

In  the  previous  cluqtter,  three  q)proximate  modds  of  the  soil  continuum  were  presented. 
The  V  «  0,  Oy  ~  0,  and  the  proposed  modds  were  rqwesented  1^  different  eqnossions  for  the 
coeffident  0  in  Eq.  2.5.  It  is  necessary  to  examine  the  accura^  of  the  solutions  given  by  each 
modd.  Wood’s  rigorous  solution  (Wood,  1973)  will  be  used  as  the  measure  of  the  accuracy  of 
the  iq>proximate  solutions.  Wood’s  solution  is  strictly  valid  for  an  dastic  layer  of  finite  length 
retained  by  a  wall  with  a  smooth  interfiice  at  each  end.  Note  that  the  "stat^  1-g  scrfution  is  the 
solution  for  very  low  fi«quendes. 

From  Eq.  3.25  the  static  deflection  produced  by  a  1-g  static  force  may  be  obtained  by 
letting  the  exciting  frequency  o  approach  zero. 


(3.28) 


The  corre^nding  1-g  static  thnist  is  obtained  by  substitutii^  Eq.  3.28  into  Eq.  3.23 
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Pm  =  pPgZZ 


b.(bj+0ai,) 


Pm  =  PPgSZ 


_ 16 _ 

jt(2n-l)*oiB  L/H 


The  1-g  static  moment  acting  at  the  base  of  the  waU  is 


Mrt  =  ppgZZ 


amamn-««(bnH) 

bi(b5  +eai) 


(3.29a) 


(3.29b) 


(3.30) 


The  total  thrust  against  the  wall  due  to  1-g  static  horizontal  force  is  determined  by  doing  a  double 
summation  for  modes  m  and  n  in  Eq.  3.29.  A  normalized  thrust  ratio  is  defined  as 


Thrust  Ratios 


Total  Thrust 
pH^  A,^ 


(3.31) 


\^ere  Am  is  the  peak  ground  acceleration  in  m/sec^  ft/sec^  or  other  consistent  units. 

The  three  ^>proxiniate  modds  are  used  to  obUun  the  total  1-g  static  force  for  two 
different  wall-soil  systems  with  L/H  *  S.O  and  L/H  *■  l.S. 
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Discnssion  of  Results 

The  results  of  the  analyses  are  compared  with  Wood’s  (1973)  solutions  for  L/H  *  5.0  and 
LTH  »  1.5  are  shown  in  Figs.  3.2(a)  and  3.2(b),  re^)ectivdy.  The  following  observations  may  be 
nuule 


•  For  both  L/H  -  5.0  and  L/H  1.5  for  ail  values  of  v,  the  proposed  modd  gives 
results  that  are  in  very  good  agreement  with  the  exact  solutions.  The  modd 
works  better  for  walls  retaining  finite  backfill  (L/H  =  1.5).  The  modd  gives  a 
total  force  slightly  less  than  the  exact  force. 

•  For  L/H  -  5.0,  the  Oy  =  0  modd  yidds  results  that  are  in  very  good  agreement 
with  the  exact  solution  for  all  v.  For  walls  with  long  backfills,  the  accuracy  of 
the  Oy  =  0  modd  is  comparable  to  that  of  the  proposed  modd.  The  Oy  =  0 
model  overestimates  the  response  by  about  8%,  the  proposed  modd 
underestimates  the  re^nse  by  about  5%. 

•  For  L/H  »  1.5,  solutions  by  the  proposed  model  are  much  closo*  to  the  exact 
solution  for  all  values  of  v.  The  Oy  -  0  model  overestimates  the  total  force  by 
up  to  18%.  The  proposed  model  underestimates  the  total  force  by  less  than 
4%. 


The  V  =  0  modd  gives  good  accura^,  provided  values  of  v  <  0.3.  As  v 
exceeds  0.3,  the  solutions  start  to  deviate  fi'om  the  exact  solutions.  For  L/H  » 
5.0,  the  accuracy  of  the  solution  becomes  unacceptable  as  v  >  0.4. 


THRUST  RATIO  THRUST  RATIO 
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fig.  i.r.  Comparison  of  approximate  solutions  for  ripd-wall  systems  with  Wood’s  (1973)  exact 

solution  for  (a)  L/H-5.0  and  (b)  L/H-1.5. 


21 


Recently,  Veletsos  and  Younan  (1994)  also  evaluated  the  accuracy  of  the  Oy  *  0  and  v  *  0 
models  compared  to  Wood’s  rigorous  solution.  Thdr  findings  regarding  these  models  agree  with 
those  above.  However,  these  more  recent  studies  were  limited  to  semi-infinhe  backfills. 

The  studies  indicate  that  the  proposed  model  gives  the  best  iq)proximation  of  solutions  for 
the  rigid-wall  ^sterns  with  elastic  backfills.  Therefore,  this  proposed  model  will  be  used  for  all 
further  studies. 
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CHAPTER4 


ELASTIC  DYNAMIC  RESPONSE  OF  PILES 


Equations  of  Motion  for  Soil-Pile  System 

The  simplified  modd  of  the  soil  continuum  derived  earlier  is  now  extended  to  incmporate 
the  important  aspects  of  3-D  response  and  adapted  to  accommodate  piles. 

Under  vertically  propagating  shear  waves  the  soils  mainly  undergo  shear  deformations  in 
horizontal  planes  except  in  the  area  near  the  pile  v^ere  extensive  compressive  defomuUions 
develop  in  the  direction  of  shaking.  The  compressive  deformations  also  generate  additional 
shearing  deformations  along  the  sides  of  the  pile  because  of  the  limited  extent  of  the  pile  cross- 
section.  In  light  of  these  observations,  assumptions  are  made  that  dynamic  motions  of  soils  are 
governed  by  two  kinds  of  shear  stresses,  and  compresdve  stresses  in  the  direction  of  shaking,  y. 
Let  v  represent  the  horizontal  displacement  of  the  soil  in  the  direction  of  shaking,  y  (Fig. 

4.1).  The  inertial  force  is  p, — The  compressive  force  in  the  direction  of  shaking  is  OG — 7^. 

St”*  ay-* 

a^v 

The  shear  waves  propagate  in  the  z  direction.  The  shear  force  in  the  xOy  plane  is  G — 7  >  ^  ^ 

dxr 


shear  force  in  yOz  plane  is  G — 7 .  .^plying  dynamic  equilibrium  in  the  y-direction,  the  dynamic 

dxr 

governing  equation  under  fi-ee  vibration  of  the  soil  continuum  is  written  as 


a^  ^  a^v  a^v 
ax^  d^  ’  at^ 


(4.1) 
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Fig.  4. 1 .  Important  aspects  of  pile-soil  interaction  during  horizontal  exdtation  by  diear  waves. 
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where  6  2/(l-v),  G  is  the  complex  shear  modulus,  and  p.  is  the  mass  density.  Since  soil  is  a 

hystoxtic  material,  the  complex  shear  modulus  G  is  expressed  as  G  »  Go  (l-H  2X),  in  which  Go  is 
the  shear  modulus  of  soil,  and  X  is  the  equivalent  viscous  damping  ratio.  Radiation  damping  will 
be  considered  later. 

Piles  are  modelled  using  ordinary  Eulerian  beam  theory.  Only  the  bending  momoit  in  the 
plane  of  shaking  (yOz)  is  included.  This  consideration  is  appropriate  in  the  case  of  earthquake 
loading  for  single  piles. 

The  dynamic  equation  of  motion  for  a  vertical  pile  is  written  as 


EpI 


a^v  d\ 


(4.2) 


where  Epiy  is  the  flexural  rigidity  of  the  pile  in  the  direction  of  shaking,  and  Pp  is  the  mass  density 
of  the  pile. 


Finite  Element  Formulation  of  Equations  of  Motion 

The  3-dimensional  soil  continuum  is  divided  into  a  number  of  3-D  finite  elements  of  the 
type  shown  in  Fig.  4.2.  The  displacement  field  at  any  point  in  each  element  is  specified  by  the 
nodal  displacements  and  appropriate  shape  functions.  A  linear  displacement  field  in  the  soil 
element  is  used  because  of  its  simplicity  and  effectiveness. 

The  displacement  vector  v(x,y,z)  is  expressed  as 
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Z 


SOIL  ELEMENT 


Z 

□  V2,(dv/dz)2 


1#1  VI.  (dv/d2)l 


BEAM  ELEMENT 


Fig.  4.2:  Finite  elements  used  for  soil  continuum  and  pile. 


v(x,y,z)  =  ZNi  X  Vi,...  i «  1,8  (4.3) 

wdiere  Ni  is  a  shape  fiinction,  and  Vi  is  a  nodal  displacement. 

Applymg  the  Galerkin's  wdghted  readual  procedure  to  Eq.  4.1,  the  complcK  valued 
stiffiiess  matrix  of  the  soil  element  is  determined 
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dtij 

dy  dy 


+  G 


aNj  gNj 

dx.  dx 


+  G 


aNj 

dz  dz  J 


dxdydz 


(4.4) 


For  a  vertical  beam  dement,  (Fig.  4.2)  the  stififoess  matrix  is 


P^lpile  ~ 


12 

6t 

-12 

6t 

Epiy 

6/ 

4/^ 

-6£ 

-12 

-6t 

12 

-6t 

6/ 

-6i 

4t^ 

(4.5) 


The  two  nodes  in  a  vertical  beam  element  are  shared  by  adjacent  soil  dements,  thus  coupling  the 
pile  and  the  soil.  The  stiffiiess  of  a  node,  therefore,  is  comprised  of  stiffiiess  contributions  from 
both  the  pile  and  the  soil. 

Consistent  mass  matrices  are  used  for  both  soil  dements  and  pile  elements.  The  consistoit 
mass  matrix  for  a  soil  element  is 


[MJsoa  =  Ps  -NiNj  - dxdydz 


The  consistent  mass  matrix  for  a  pile  element  is  written  as 


(4.6) 


IMjpfle  = 


156 

22/ 

54 

-13/ 

PpA/ 

22/ 

4/2 

13/ 

-3/2 

420 

54 

13/ 

156 

-22/ 

-13/ 

-3/2 

-22/ 

4/2 

(4.7) 
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The  global  stifBiess  matrix  [K]  and  the  ^obal  mass  matrix  [M]  are  constructed  by  standard 
routines. 

The  radiation  damping  is  modelled  by  using  a  set  of  set  of  dashpots  along  the  pile  shaft. 
The  damping  force  F4  per  unit  length  along  the  {Hie  is  considered  proportional  to  the  vdodty  and 
is  given  by 


Fd 


(4.8) 


where 


c,  =  6p.V.<la5“'‘* 


(4.9) 


where  c,  is  the  radiation  dashpot  coefficient. 

Gazetas  et  al.  (1993)  proposed  simple  e}q)resdons  for  the  radiation  dashpot  coefficients  c, 
for  both  horizontal  motion  and  vertical  motion.  The  element  radiation  damping  matrix  [C]«h.  is 


[^lelem  “  ^20 


156 

22t 

54 

-13/ 

221 

13/ 

-3/2 

54 

13/ 

156 

-22/ 

-m 

-3/2 

-22/ 

4/2 

(4.10) 


The  global  dynamic  equilibrium  equation  in  matrix  form  is  now 


[M]{v}+[C]{v}+[K]{v}  =  {P(t)} 


(4.11) 
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in  which  {v},  {v}  and  {v}  are  the  nodal  accelerations,  vdodties  and  diq>lacement8^  re^Mctivdy, 
{P(t)}  ve  external  dynamic  loads. 

Solutions  of  Eq.  4.10  will  now  be  devdoped  to  ^ve  the  stiffoess  and  daniq>ing  oi  single 
piles  (pOe  impedances)  as  a  function  of  frequency.  These  impedances  will  be  conq>ared  with 
impedances  calculated  using  the  full  3-D  equations  for  the  continuum  to  verify  the  accuracy  of  the 
proposed  modd  for  pile  analysis. 
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CHAPTERS 

PILE  IMPEDANCES:  SOLUTIONS  FOR  HARMONIC  LOADING 

The  impedance  Ky  are  defined  as  the  amplitudes  of  harmonic  forces  (or  moments)  that 
have  to  be  applied  at  the  pile  head  in  order  to  generate  a  harmoiuc  motion  with  a  unit  amplitude  in 
the  direction  of  the  specified  degree  of  fi’eedom  (Novak,  1991)  as  shown  in  Fig.  5.1. 


Fig.  S.l:  Pile  head  impedances. 

The  impedances  are  defined  as  follows: 

Kw:  the  complex-valued  pile  head  horizontal  force  required  to  generate  unit  bHrimitfui 
displacement  (u-1.0)  at  the  pile  head  while  the  pile  head  rotation  is  fixed  (6^). 
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K««;  the  complex-valued  pile  head  moment  generated  by  the  unit  lateral  diqdaoement  (u*1.0) 
at  the  pile  head  whQe  the  pile  head  rotatimi  is  fixed  (0-0) 

Km:  the  complex-valued  pile  head  moment  required  to  generate  the  unit  pile  head  rotation 
(O-'l.O)  while  pUe  head  lateral  di^lacement  is  fixed  (u-O) 

Since  the  pfle  head  impedances  K^,  Km,  Km  are  conq)lex  valued,  they  are  usually 
expressed  by  their  real  and  imaginaiy  parts  as 

Kjj  =k,3 +iC5  =k^ +i(DCjj  (5.1) 

in  which  kg  and  Cg  are  the  real  and  imaginaiy  parts  of  the  complex  impedance,  reqiectivdy,  and 
are  usually  referred  as  the  stifiBiess  and  damping  at  the  pile  head,  i  s  c^  =  /o>  > 

coefiBciem  of  equivalent  viscous  damping;  and  co  is  the  drcular  firequency  of  the  iqiplied  load.  All 
the  parameters  in  Eq.  S.l  are  dependent  on  fi-equency  q>.  Determination  of  the  pile  impedances 
requires  solutions  of  the  equations  of  motion  for  harmonic  loading. 

Under  harmonic  loading  P(t)  *  Poe***,  the  displaconent  vector  is  of  the  form  v  *  voe***,  and 
Eq.  4.1 1  is  rewritten  as 


{[Kl  +  i.a)[C]-a)2[M]){Vo}  =  {Pg} 


(5.2) 


or 


[Klgiobd  (Vo)  =  (Po) 


(5.3) 
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where 


fKJgioW  »  rK]+i  .«(C]-®*CMl 


(5.4) 


The  |Hle  head  tiiq>edances  IC««  and  Km  can  be  obtained  using  Eq.  4.8 ^>ldying  ^>propriate 
loading  and  fixity  conditions  at  the  pfle  head. 


Determination  of  Impedances  Kr,  and 

If  the  pile  head  is  fixed  against  rotation  and  a  unit  horizontal  di^lacement  ai^lied  at  the 
pUe  head,  then  Eq.  5.3  has  the  form 


La] 

Vq 

*  ^ 

0 

[^Iglofaal 

1.0 

*  ^  * 

Kw 

0.0 

c  J 

.K«0. 

(5.5) 


vAten  are  the  di^lacements  of  the  nodes  other  than  pile  head.  Eq.  5.5  can  be  solved  by 
dividing  the  equation  by  K«v  and  elimirudiiig  tlw  row  correspondirig  to  zero  rotation 


K]g|otMl 


(5.6) 


Tlus  equation  suggests  that  an  easier  ahemative  for  detemunirig  inqredances  is  to  ^ipiy  a  unit 
horizontal  fierce  at  the  pile  head  and  calculate  the  complex  displacement  vo^  The  moment  at  pie 
head  Mo'  corre^nding  to  vo'  is  also  calculated. 
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Then  the  pUe  head  in^edance  Kw  is  deteraune^ 


K 


YV  ~ 


(5.7) 


and  the  correqranding  cross>mq)edance 


(5.8) 


Determination  of  Rotational  Impedance  Km 

For  a  pile  head  fixed  against  horizontal  displacement  and  given  a  unit  pile  head  rotation 
Eq.  5.3  has  the  form 


V? 

0 

^Iglobd 

1.0 

»  ^  < 

Kev 

.00 

.Kee. 

(5.9) 


(5.10) 


vdiere  vo^  are  the  displacements  of  nodes  other  than  pile  head.  The  pfle  head  rotation  Oo'  is  tiie 
complex  response  to  the  application  of  a  unit  moment  at  the  pile  head. 
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The  rotatioiial  iii^>edance  K«  is  then  9ven  as 

Ko9  -  ^  (511) 

e; 


Because  of  redprodty  prindple,  the  cross-coupling  impedances  and  Ka.  are  identical. 

Dynamic  Impedance  of  Single  Piles 

Anal^es  are  conducted  to  determine  the  frequency  dependence  of  the  inq>edances  K«v, 
Kve  and  Km  for  single  piles.  The  dimensionless  frequency  ratio  ao,  is  used  to  characterize  the 
frequency,  vdiere  ao  is  defined  by 


in  which  co  is  the  angular  frequency  of  the  exdting  loads  (force  or  moment)  at  the  pile  head,  d  is 
the  diameter  of  the  pile,  and  V,  is  the  shear  wave  velodty  of  the  soU  medium.  For  a  uniform  soil 
profile  with  a  shear  modulus  G  and  a  mass  denaty  p,  V«  »<^G7p. 

For  given  values  of  Ep/E,  and  ao,  the  ratios  Kv«/(E|d),  and  Kee/(E,d^),  remain 

constam  as  the  soil  modulus  E«  varies  if  the  soil  profile  is  uniform.  Tho'efore  the  normalized 
inqmdances  KM/(E,d^,  and  Koe/(Ei(P)  are  used  in  presenting  the  fi^ueni^  dq>endaice 

of  the  impedances. 
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To  assess  the  accuracy  of  the  proposed  nmplified  3-D  finite  dement  approach,  the 
impedance  functions  of  the  pile-soil  system  shown  in  Fig.  S.2  were  computed  over  a  range  of 
fi-equencies  and  compared  with  the  full  3-D  linear  dasdc  solutions  of  Kaynia  and  Kausd  (1982). 
However,  the  latter  solutions,  though  regarded  as  benchmark  solutions,  are  not  exact.  They 
assume  the  form  of  the  pressure  distribution  between  soil  and  pile  at  any  devation  to  be  uiuform 
around  the  laterally  displacing  pile. 

single  piles: 

Ep/Es  «  1,000 
L/d  >  IS  (floating  pile) 
soil  damping  5% 
poisson's  ratio  0.4 


L 


Fig.  5.2;  A  pile-soil  system  used  for  computing  impedance  functions. 

The  impedances  are  determined  for  Ep/E. «  1,000,  Poisson’s  ratio  v  «  0.4,  and  a  damping 
ratio  X  =  5%. 

The  finite  element  mesh  shown  in  Fig.  5.3  is  used  for  the  analysis.  The  heavy  dark  verticd 
dement  is  the  pile.  Because  of  the  symmetric  character  of  the  problem,  only  half  of  the  fuU  mesh 
is  required  to  model  the  response  of  pile-soil  interaction.  The  half  mesh  consists  of  1463  nodes 
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Fig.  5.3 :  Finite  element  modelling  of  angle  pile  for  computing  impedance. 


and  1089  elements.  The  use  of  symmetry  reduces  the  size  of  the  global  matrix  by  a  factor  of  4.0 
and  greatly  reduces  the  computational  time.  The  computing  time  is  300  seconds  for  each 
frequency  using  a  33  Hz  486  PC  computer. 

The  half-mesh  includes  only  half  the  pile  so  that  for  the  linear  beam  element  Ep  and  Pp 
must  be  reduced  by  a  factor  of  2.  The  loads  applied  to  the  pile  must  also  be  reduced  from  1.0  to 
0.5.  The  effects  of  these  reductions  must  be  kept  in  mind  when  the  corresponding  impedances  are 
calculated. 

Discussion  of  Results 

Since  the  normalized  impedances  are  complex  quantities,  they  are  given  in  terms  of  thdr 
real  (stiffitess)  and  imaginary  (damping)  parts.  Tire  computed  normalized  stiffiiesses  and  damping 
ratios  as  functions  of  dimensionless  frequency  ao  are  compared  imth  the  solutions  of  Kaynia  and 
Kausel  (1982)  in  Figs.  5.4, 5.5  and  5.6. 
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In  general,  the  impedances  computed  by  the  simplified  method  agree  wdl  with  those  of 
Kaynia  and  Kausel  (1982).  The  horizontal  stifihiesses  kw  computed  by  Kaynia  and  Kausel  (1992) 
are  about  10%  larger  than  those  computed  using  the  simplified  model.  The  other  two  stifiBiesses 
kve  and  kee  show  closer  agreement.  The  differences  are  on  the  average  about  5%. 

It  is  hard  to  decide  which  of  the  two  solutions  represents  more  closely  the  true  solution  of 
this  problem  because  both  were  obtained  by  iq;)proximate  methods.  In  Kaynia  and  Kausel's 
solution  the  horizontal  load  applied  to  soil  medium  by  the  pile  is  assumed  uniformly  distributed  on 
the  cylindrical  soil  surface  around  the  pile.  This  does  not  model  well  the  non-uniform  distribution 
of  lateral  pressure  on  the  cylindrical  surface  of  the  surrounding  soil  caused  by  the  laterally 
deflecting  pile 

The  size  and  rumber  of  the  finite  elements  affect  the  computed  value  of  the  impedance. 
The  accuracy  increases  as  the  number  of  finite  elements  increases  especially  as  the  fi-equency 
increases.  At  high  firequencies  (such  as  ao  >  0.4),  large  numbers  of  finite  elements  are  needed  to 
capture  the  number  of  modes  that  are  significant  for  response  at  that  fi’equency. 

Figure  5.7  shows  a  comparison  of  the  dynamic  stiflhiesses  computed  by  two  different 
meshes.  It  is  clear  that  mesh  size  becomes  sigiuficant  when  the  fi’equency  becomes  high.  A  finer 
mesh  is  needed  to  represent  the  dynamic  responses  accurately. 

For  earthquake  loading,  the  non-dimensional  fi'equency  for  pile  foundations  is  usually  less 
than  a«  =  0.3  for  the  important  fi-equencies  in  the  ground  motions.  In  this  fi’equency  range,  the 
approximate  method  proposed  is  very  accurate,  even  with  a  rather  coarse  mesh. 
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CHAFrER6 

NONLINEAR  DYNAMIC  RESPONSE  ANALYSIS  OF  Pnj:S 

The  equations  of  motion  are  given  by  the  incremental  form  of  Eq.  4.8  as 

M]m  +[C]{Av}+[K]{v)  =  {P(t))  =  -[M][I]  AVfc(t)  (6.1) 

Analysis  of  nonlinear  response  must  be  conducted  in  the  time  domain.  The  direct  step  by  step 
integration  procedure  developed  by  >^son  et  al.  (1973)  is  used  to  integrate  the  equations  of 
motion. 

Rayldgh  damping  is  used  to  model  the  hysteretic  damping  of  the  soil  for  nonlinear 
analyds.  The  damping  element  matrix  is  ffven  by 

tcid„  =  (6.2) 

in  which  a  and  P  are  constants  related  to  the  wscous  damping  ratio  for  the  element.  Let 

a  =  ^dem  and  p  =  (6-3) 

vdiere  X«i«ni  is  the  damping  ratio  corresponding  to  element  shear  strainand  Oi  is  the  fimdamoital 
frequency  of  the  system  (Idriss  et  al.,  1974). 

The  global  damping  matrix  [C]  is  the  aggregate  of  all  the  element  damping  ratios  and  the 
radiation  damping  elements  along  the  pile. 
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For  nonlinear  analysis  because  of  integration  in  the  time  domain,  it  is  more  cfBdent  to  use 
a  diagonal  mass  nutrix  rather  than  the  consistent  nuus  matrix  used  eariier  for  harmonic  analysis. 
Th««fore,  the  mass  matrices  for  soil  and  beam  dements,  [M]Ma  and  [M]bM,  req>ectivdy,  are 
given  by 

[M]*,a  =  {10. 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0)  (6.4) 

and 

Mlbcam  =  Pp  Af  {1/2,1/78,1/2,1/78}  (6.5) 

Soil  modulus  and  damping  in  soils  are  shear  strain  dependent  (Seed  and  Idriss,  1970). 
During  analysis  compatibility  is  maintained  between  the  computed  shear  strains  and  the  effective 
modulus  and  damping  in  each  fimte  element.  The  compatibility  can  be  restored  for  each  time 
increment  during  integration  of  the  equations  of  motion  or  at  specified  times  which  are  multiples 
of  the  time  increment  for  integration.  This  procedure  differs  fi'om  the  equivalent  linear  method 
used  in  programs  such  as  SHAKE  (Schnabel  et  al.,  1972)  in  vdiich  compatibility  is  enforced  only 
after  the  complete  response  analyris  has  been  completed.  Ensuring  final  compatibility  in  that  case 
requires  iterative  anatysis  using  the  entire  duration  of  the  earthquake  in  each  aiuilysis.  No  iterative 
analyses  are  required  when  compatibility  is  enforced  during  the  analysis. 

Two  other  features  distinguish  the  nonlinear  model  proposed  here  fi’om  the  Schnabel  et  al. 
(1972)  model.  Shear  yielding  is  incorporated  by  introducing  a  very  low  modulus  \riien  the 
strength  of  the  soil  is  reached.  No  tensile  stresses  are  allowed.  This  is  accomplished  by 


introducing  a  very  low  modulus  when  the  normal  stress  in  ai^  direction  tends  to  become  greater 
than  the  tensile  strength  of  the  soil  if  any. 

The  nonlinear  method  of  analysis  will  be  vahdated  now  using  data  from  strong  shaking 


tests  of  a  ^gle  pile  in  a  centrifuge. 


44 


CHAPTER? 

CENTRIFUGE  TESTS  ON  A  SINGLE  PILE  UNDER  STRONG  SHAKING 

Test  Set-Up 

The  nonlinear  method  is  used  to  analyze  the  setsmic  response  of  a  sii^  pile  in  a 
centrifuge  test  which  was  carried  out  on  tl^  California  Institute  of  Tedinology  (Caltech) 
centrifuge  by  B.  Gohl  (1991).  Details  of  the  test  may  also  be  found  in  a  papa*  by  Rnn  and  GoM 
(1987).  Fig.  7.1  shows  the  soil-pile-structure  syston  used  in  the  test.  A  209.S  mm  long  sUunless 
steel  tube  pile  having  an  outade  diameter  of  9.52  mm  and  a  wall  thickness  of  0.25  mm  is 
embedded  in  a  dry  loose  sand  fovmdation.  The  model  pile  is  instrumented  by  8  pairs  of  foil  type 
strain  gauges  mounted  on  the  outside  of  the  pile  to  measure  bending  strains  at  the  locations 
shown  in  Fig.  7. 1.  An  average  centrifuge  accel«ation  of  60g  was  used  in  the  t^. 

The  pile  has  a  free  standing  length  of  16.5  mm  above  the  soil  sur&ce.  The  effect  of  a 
super-structure  is  simulated  by  clamping  a  rigid  mass  to  the  head  of  the  pile.  The  wdght  of  the 
structural  mass  including  the  pile  head  insert  and  the  pile  head  clamp  is  2.416  N.  The  noass 
moment  of  inertia  about  the  centre  of  gravity  is  I,,  =  0.0683  N.sec^.nun.  The  centre  of  gravity  of 
the  mass  is  located  16.5  mm  above  the  pile  head.  The  model  pile  has  an  average  flexural  rigidity 
of  13.26  N.m^  and  a  mass  density  of  74.7  kN/m^. 

The  pile  head  mass  is  instrumented  uai^  a  non-contact  photovoltaic  di^lacement 
transducer  and  an  Entran  miniature  accelerometer.  The  locations  of  the  accelerometer  and  light 
emitting  diode  (L.E.D.)  used  by  the  displacement  sensor  are  shown  in  Fig.  7.1.  The  pile  head 
displacements  are  measured  with  respect  to  the  moving  base  of  the  soil  container. 
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Fig.  7.1.  The  layout  of  the  centrifijge  test  for  a  single  pile. 

The  sand  used  for  the  test  was  a  loose  sand  with  a  void  ratio  e»  *  0.78  and  a  mass  density 
p  «  1.50  Mg/m^  Gohl  (1991)  has  shown  that  the  low  strain  shear  moduli  of  the  sand  foundation 
vary  as  the  square  root  of  the  depth,  and  th^  can  be  quantitativdy  evaluated  using  the  Hardin  and 
Black  (1968)  Eq.  (7.1) 


I 

I 

I 


(2.973  -  e- )*  ,^  ,0^ 
G....  =  3230^  (rf,)" 


(7.1) 
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where  e*  is  the  void  ratio  of  the  sand  and  a*,  is  the  mean  nwmal  effective  confiniiy 
pressure  in  kPa. 

A  horizontal  acceleration  is  b^ut  at  die  base  of  the  system.  The  peak  acoderatioa  of  die 
input  motion  isO.lSSg.  The  computed  Fourier  ain|ditude  ratios  ofdie  pile  head  response  and  die 
free  fidd  motion  with  respect  to  the  ir^Mit  modoiis  are  given  in  Fig.  7.2(a)  and  ^  7.2(b).  The 
natural  frequency  of  the  free  fidd  acceleration  is  estimated  to  be  2.75  Hz  fiom  Rg.  7.2(b)  and  die 


Fig.  7.2.  The  Fourier  Spectra  of  acederations  (after  Gohl,  1991). 
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fimdamattal  firequency  of  the  pile  (fixMn  Fig.  7.2(e))  to  be  1.1  Hz.  The  period  ot  peak  pile 
response  is  much  Icmger  than  the  period  of  the  free  field. 

The  centrifuge  test  was  analyzed  1^  the  sinq)lified  3-D  finite  dement  method  analysis. 
Fig.  7.3  shows  the  finite  dement  modd  used  fr>r  analysis.  The  sand  dq)Osit  is  dhnded  into  11 
layers.  Li^  thickness  is  reduced  as  the  sdl  suifiice  is  approached  to  aflow  more  detdled 
meddling  of  the  stress  and  strain  fidd  vdiere  lateral  soil-|»le  interaedem  is  stremgest.  The  |Mle  is 
modelled  using  IS  beam  dements  induding  5  dements  above  the  sofl  surftoe.  The  super¬ 
structure  n»ss  is  treated  as  a  rigid  body. 


GEO. SCALE  0 - - - ^ 


Fig.  7.3.  The  finite  dement  modelling  of  centrifiige  test. 

The  filnhe  dement  analysis  is  carried  out  in  the  time  domain.  Nonlinear  analysis  is 
performed  to  accoum  for  the  changes  in  shear  moduli  and  damping  ratios  due  to  dynamic  shear 
strains.  The  shear-strain  dependency  of  both  the  shear  modulus  and  damping  ratio  used  in  the 
analysis  is  shown  in  Fig.  7.4.  The  low  strain  shear  moduli  Gm*  were  determined  using  Eq.  (7.1). 
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SHEAR  STRAIN  («) 

Fig.  7.4.  The  relationships  between  shear  modulus^  damping  and  the  shear  strain. 


Discussion  of  Results 

The  computed  acceleration  response  at  the  pile  head  is  plotted  against  the  measured 
response  in  Fig.  7.5.  Fairly  good  agreement  between  the  measured  and  the  computed 
accelerations  is  observed  in  the  r^on  of  strong  shaking. 

The  computed  time  histories  of  moments  in  the  pile  at  the  ground  sur&ce  and  at  a  dqrth  of 
3  m  are  plotted  against  the  corresponding  recorded  time  histories  in  Hg.  7.6  and  Fig.  7.7.  There 
is  satis&ctoiy  agreement  between  the  computed  and  measured  moments  in  the  range  of  laiga* 
moments.  The  computed  and  measured  moment  distributions  along  the  pile  at  the  moment  of 
peak  pile  head  deflection  are  shown  in  I^g.  7.8.  The  computed  moments  agree  quite  wdl  with  the 
measured  moments.  The  moments  increase  to  a  maximum  value  at  a  depth  of  3.5  diameters,  and 


ACCELERATI 


BENDING  MOMENT  (KN.M) 
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Fig.  7.7.  The  computed  versus  measured  moment  response  at  dq>th  D  «  3  m. 


-400  100 

BENDING  MOMENT  (KN.M) 


Fig.  7.8.  The  computed  versus  measured  moment  distribution  of  the  pile  at  peak  pile  deflection. 
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then  decrease  to  zero  at  a  dq;)th  around  12.S  diametoa.  The  momoits  along  the  pile  have  the 
same  agns  at  ai^  instant  of  time,  suggesting  that  the  inertial  intoaction  caused  by  the  inle  head 
mass  dominates  response,  and  the  pile  is  vibrating  in  its  first  mode.  The  peak  monMnt  predicted 
by  the  ^plified  3-D  finite  dement  analysis  is  344  IcNm  compared  with  a  measured  peak  value  of 
325  kNm  an  error  of  only  6%. 

It  is  interesting  to  show  the  degradation  in  shear  modulus  with  shear  strain  around  the  pile 
during  shaking.  Distributions  of  moduli  at  spedfic  depths  and  a  q)edfic  tioM  during  the 
earthquake  are  shown  in  Fig.  7.9.  The  figures  show  that  dgnificant  modulus  degradation  occurs 
near  the  pile  and  is  most  pronounced  near  the  pile  head. 

Computational  Time 

Using  a  PC-486  (33  MHz)  computer,  the  average  CPU  time  needed  to  complete  one 
integration  step  is  7  sec  for  the  finite  element  grid  diown  in  Fig.  7.3,  and  3  hours  of  CPU  time  are 
required  for  the  fiiU  input  acceleration  record  of  1550  steps.  The  computational  time  would  be 
shorter  for  a  linear  elastic  analysis. 

Pile  Impedances  During  Strong  Shaking 

Dynamic  impedances  as  functions  of  time  were  computed  corresponding  to  the  strain 
dependent  shear  moduli  fi'om  the  finite  element  aiudytis.  Harmonic  loads  with  an  amplitude  of 
unity  were  applied  at  the  pile  head,  and  the  resulting  equations  were  solved  to  obtain  the  complex 
valued  pile  impedances.  The  pile  impedances  were  evaluated  at  the  ground  surfiu:e.  This  is  the 
first  time  that  the  time  histories  of  pile  impedance  during  an  earthquake  have  been  determined. 


S2 


single  pile  at  17. 11  sec 


initial  shear  modulus  1 2945  kPa 


(a)  at  depth  0.25  m 


initial  shear  modulus  36610  kPa 


(b)  at  depth  2.10  m 


Fig.  7.9.  3-D  plots  of  the  distribution  of  efifective  shear  moduli  with  dqrth  around  a  jrile  during 

dynanuc  excitation. 
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The  dynamic  stiffiiess  of  the  pile  decreases  dramatically  as  the  level  of  shaking  inneases 
(Fig.  7.10).  The  dynamic  stiffiiesses  experience  thdr  lowest  values  between  about  10  and  14 
seconds,  when  the  maximum  accelerations  occurs  at  the  pile  head.  It  can  be  seen  that  the  lateral 
stiffiiess  component  kw  decreased  more  than  the  rotational  stiffiiess  kee  or  the  coupled  lateral- 
rotational  stiffiiess  k«e-  On  the  other  hand  the  equivalent  damping  coefficients  increased  as  the 
level  of  shaking  increased  because  the  hysteretic  damping  of  the  soil  increased  with  increased 
strains. 


Fig.  7. 10.  Impedances  kyu,  k«B  and  kee  of  the  single  pile. 


At  its  lowest  level,  kw  decreased  to  22,000  Id^/m,  only  15.2  %  of  its  initial  stiffiiess  of 
145,000  kN/m.  kve  decreased  to  45,000  kN/rad  or  36%  of  its  initial  stiffiiess  of  125,000  kN/rad. 
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kae  showed  the  least  effect  of  shear  strain.  It  decreased  to  138,000  IcN.m/rad  or  63.6%  of  its 
initial  stiffiiess  of  217,000  kN.m/rad.  The  stifihiesses  rounded  when  the  level  of  shakif^ 
decreased  with  time.  Representative  values  of  the  pile  stiffiiesses  Icw,  lc««  and  Icm  that  might  be 
used  in  structural  analyses  may  be  taken  as  40,000  kN/m,  65,000  kN/rad  and  160,000  kN.m/rad, 
reflectively.  These  stiffiiesses  are  32%,  52%  and  73.7%  of  the  original  stiffiiesses. 

The  effect  of  frequency  on  both  sdffiiess  and  damping  is  explored  for  a  range  of 
frequencies  from  1.91  Hz  to  10  Hz  at  different  times  during  the  dynamic  shaking  of  the  pile.  The 
stiffiiess  response  is  shown  in  Fig.  7. 1 1 .  Clearly  within  this  range  of  frequences  which  is  typical 
of  the  frequencies  of  peak  response  in  many  near  and  medium  field  earthquake  motions,  there  is 
no  frequency  effect  on  stiffiiess. 


Fig.  7. 1 1 .  Variation  of  pile  head  stiffiiesses  with  frequency. 
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The  damping  response  is  shown  in  Fig.  7.12.  Clearfy,  even  within  this  rdatively  small 
frequency  range  there  is  a  agnificant  variation  in  the  damping  of  the  pile.  This  indicates  that  it  is 


very  difBcuh  to  select  the  proper  equivalent  dashpot  to  reflect  the  damping  of  a  pfle  foundation  in 


a  structural  analysis. 


Fig.  7. 12.  The  effect  of  frequency  on  pile  damping. 


The  centrifuge  test  provides  an  opportunity  to  evaluate  the  effects  of  inertial  interaction  on 
the  stiffiiess  of  a  pile  foundation.  Many  procedures  in  practice  for  evaluating  pile  stiffiresses  are 
based  on  computing  the  inertial  and  kinematic  components  of  soil  structure  interaction  separately 
(Gazetas,  1991a,  1991b).  This  is  acceptable  for  elastic  response  because  the  additionai 
foundation  excitation  caused  by  the  inertia  of  the  structural  mass  does  not  affect  the  stifl&iess  of 
the  foundation.  However,  when  there  is  nonlinear  response  of  the  foundation,  the  inertial 
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interaction  of  the  structure  with  the  foundation  can  cause  major  changes  in  foundation  stiflhess 
and  in  the  period  of  peak  response. 

This  is  clearly  shown  by  analytis  of  the  angle  pile  in  the  centrifoge  test  with  and  without 
the  structural  mass.  The  results  of  the  analysis  are  shown  in  Fig.  7. 13  '^ch  shows  the  significant 
degradation  in  pile  stiffiiess  due  to  inertial  mass  at  the  pile  head  for  both  translational  and 
rotational  stifiGhesses,  kw  and  respectively.  Cleaiiy,  for  strong  earthquake  shakitig  the  effects 
of  inertial  nuiss  carmot  be  ignored  and  relying  on  kinematic  stififoess  only  may  lead  to  a  serious 
overestimation  of  pile  stifihiess. 

So  for,  validation  of  the  simplified  method  of  analysis  has  been  done  by  comparing 
solutions  with  published  elastic  solutions  using  full  3-D  formulations  or  by  data  fi'om  centrifoge 
tests.  In  conclusion,  a  forced  vibration  test  on  a  foil  size  pile  in  the  field  will  now  be  atudyzed. 


STIFFNESS  STIFFNESS 


translational  sH/fnoss  kyj(kN/m) 


Fig.  7.13.  The  effect  of  inertial  interaction  on  pile  head  stiffiiesses. 
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CHAPTERS 

ANALYSIS  OF  FULL-SCALE  VIBRATION  TEST  ON  A  FRANK!  PILE 


Test  Set-Up 

A  full  scale  vibration  test  on  a  single  Frank!  type  pile  was  performed  by  Sy  and  Shi  (1992) 
at  the  University  of  British  Columbia  Pile  Research  Facility  located  in  the  Fraser  river  ddta  just 
south  of  Vancouver,  B.C..  The  soil  profile  at  the  testing  site  consists  of  4  m  of  sand  and  gravd 
fill  overlying  aim  thick  silt  layer  over  fine  grained  sand  to  a  depth  of  40  m.  A  seismic  cone 
penetration  test  (SCPT  88-6)  was  conducted  0.9  m  from  the  test  pile.  In  addition,  SPT  tests  were 
conducted  in  a  mud-rotary  drill  hole  (  DH88-2)  2.4  m  from  the  test  pile.  The  measured  in-situ 
shear  wave  velocity  data  are  presented  in  Fig.  8.1,  together  with  the  cone  penetration  test  (CPT) 
data  and  the  Standard  Penetration  Test  (SPT)  data. 

The  layout  of  the  pile  test  is  shown  in  Fig.  8.2.  The  pile  has  an  oqianded  ^herical  base 
with  a  nominal  diameter  of  0.93  m.  For  6.4  m  above  the  expanded  base,  the  pile  has  a  nomiiud 
diameter  of  510  mm.  The  remaining  length  has  a  square  cross-section  with  a  side  width  of  510 
mm.  A  structural  mass  consisting  of  1.6  m  cube  of  rdnforced  concrete  was  formed  on  top  of  the 
pile  with  a  clearance  of  150  nun  above  the  ground  surfru^e. 

Accelerometers  were  mounted  on  the  shake  mass  and  the  pile  cap  to  measure  the  (fyruunic 
input  force  and  the  pile  cq>  responses  as  shown  in  Fig.  8.2.  The  vertical  and  coupled  horizontal 
and  rocking  modes  of  vibration  were  obtained  by  rotating  the  shaker  so  that  the  dyiuunic  fisrces 
were  applied  in  the  vertical  and  horizontal  directions.  The  natural  fiaquency  of  the  aq)-pfle-soil 
system  in  each  vibration  mode  was  estimated  by  applying  random  bandwidth  excitation.  Then  a 


Fig.  8.1.  The  in-^  measured  geotechnical  data. 


Fig.  8.2.  The  layout  of  the  full-scale  vibration  test  on  a  angle  pile. 
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detailed  sinuscHdal  frequency  iweep  was  carried  out  around  the  natural  frequency  indicated  frtmi 
the  random  bandundth  test. 

The  resonant  frequoicies  from  the  anusoidal  sweep  testing  are  evident  in  Fig.  8.3.  The 
damping  ratios  are  calculated  fit)m  the  measured  frequency  re^nse  curves  shown  in  Fig.  8.3 
uang  the  half  power  point  or  bandwidth  method  (Clough  and  Penzien,  1975). 


nteouENor  (hs) 


FREQUENCY  (Hs) 


Fig.  8.3.  The  measured  dynanuc  responses  of  the  {rile  cap  under  rinusoidal  iiqnit 

(after  Sy  and  Siu,  1992). 
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Analysis  of  Teits 

The  structural  properties  of  the  pile  cap  and  the  i  {ale  used  in  the  analysis  are  presented 
in  Table  8.1.  The  slwar  wave  vdodty  V.,  unit  weight  y,  and  daiiq>ing  ratio  used  in  the  anahrsis 
are  shown  in  IHg.  8.4.  Follov^  Sy  and  Siu  (1992X  an  upper  bound  of  the  measured  V,  values 
was  used  to  account  for  the  effect  of  soil  densification  caused  by  pile  installation  excqrt  fiar  the 
top  1.2  m  depth.  V«  values  in  the  upper  1.2  m  were  reduced,  since  the  (viginal  soH  around  the 
extended  pUe  shaft  section  was  replaced  by  the  loose  backfill.  A  Poisson's  ratio  of  0.3  is  assumed 
for  all  soil  layers. 


Table  8.1.  Structural  properties  of  pile  cap  and  test  pile  (after  Sy  and  Siu,  1992). 


PARAMETER 

UNIT 

VALUE 

PILE  CAPE  AND  SHAKER 

Mass 

Mg 

10.118 

Mass  moment  of  inertia 

Mgm^ 

4.317 

Hdght  to  centre  of  gravity 

m 

0.8 

TEST  PILE 

Top  1.37m:  axial  rigidity  (EA) 

MN 

6350 

Top  1 .37m:  flexural  ri^dity  (El) 

141 

1.37-7.77m:  axial  rigidity  (EA) 

MN 

5150 

1.37-7.77m:  flexural  rigidity  (El) 

92 

Base:  axial  rigidity  (EA) 

MN 

14,720 

Base:  flexural  rigidity  (EA) 

MNm’ 

800 

Material  damping  ratio 

0.01 

Poisson’s  ratio 

0.25 

62 


UNir  wr 
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Fig.  8.4.  The  soil  parameters  used  in  analysis  (Sy  and  Siu,  1992). 

Figure  8.S  shows  the  3*D  finite  element  model  used  in  present  analyas.  Due  to  symmetry, 
only  half  of  the  full  mesh  is  needed.  The  finite  dement  modd  condsts  of  1225  nodes  and  889 
dements.  There  is  one  beam  element  above  the  ground  surface  to  represem  the  pfle  s^ment 
above  the  ground.  The  expanded  base  is  modelled  by  a  solid  element  rather  than  a  beam  dement 
in  the  finite  element  analysis. 

Since  the  pile  behaves  elasticdly  under  the  veiy  low  excitation  forces  used  in  the  test  (Fm 
*  165  N),  it  is  possible  to  use  an  uncoupled  analyds  and  treat  the  pile  foundation  and  the  pUe 
structure  above  the  ground  separately.  First  the  pile  impedances  are  determined  as  a  fiinction  of 
frequency.  Then  the  response  of  the  mass  on  the  pile  head  is  computed  incorporating  the  proper 
stiffiiess  and  damping  components  of  the  impedances  depending  on  the  fi^uency  of  the  exciting 


force. 
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Fig.  8.S.  Finite  element  model  used  in  analysis. 


The  pile  head  impedances  (stifi&iess  and  damping)  of  the  pile  foundation  were  con^Hited 
using  the  proposed  simplified  3-D  finite  element  method.  Harmonic  force  or  moment  with  unit 
amplitude  was  applied  at  the  pile  head,  and  the  resulting  complex  valued  displacement  at  the  pile 
head  was  determined.  Impedances  were  then  evaluated. 

The  horizontal  and  rocking  responses  of  the  pile  head  nwss  were  obtained  by  using  the 
two-degree  of  fi-eedom  system  shown  in  Fig.  8.6(b).  The  coupled  translation-rotation  equation  of 
motion  in  £q.  8.1  describes  the  motion  of  the  ^em. 


m  m'hpg  Vp  ^  kyy +iCyy  ky9+iCyg  f'^pl_/^o 
m*hcg  leg  .  .®p,  .^ty9+iCye  kgg+iCge.^OpJ  IMq. 


where  m  is  the  nnass  of  the  pile  cap  and  shaker,  heg  is  the  height  of  the  centre  of  gravity  to  the  pile 
head;  leg  is  the  mass  moment  of  inertia  at  the  centre  of  gravity;  ky  and  C,  are  the  stifi&iesses  and 
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Fig.  8.6.  Uncoupled  systems  for  modelling  the  vertical  motions  (a),  and  the  translational  and 

rotational  motions  (b). 


damping  at  the  pile  head;  Vp  and  Op  are  the  translation  and  rotation  at  the  pile  head;  and  Po  and  Mo 
are  amplitudes  of  the  harmonic  external  force  and  moment,  respectively,  applied  at  the  pile  head. 

Harmonic  load  amplitudes  of  165  N  in  the  horizontal  direction  with  a  coupled  moment  of 
335  N.m  were  applied  at  the  pile  head  to  simulate  conditions  created  by  the  test. 

The  horizontal  displacement  amplitude  at  the  centre  of  gravity,  Ve,  is  given  by 


Vcg=Vp+ehcg  (8.2) 

The  analyses  were  carried  out  at  different  frequencies  o.  The  computed  horizontal 
displacement  amplitude  versus  fi'equency  (o  is  shown  in  Fig.  8.7.  A  very  dear  and  pronounced 
peak  response  is  observed  at  a  fr-equency  of  6.67  Hz  compared  to  the  measured  frequency  of  6.5 
Hz.  The  computed  and  measured  frmdamental  firequendes  and  damping  ratios  for  the  trandation 
and  vertical  modes  of  response  are  9ven  in  Table  8.2.  The  agreement  between  them  is  very  good. 
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Fig.  8.7.  Amplitudes  of  horizontal  displacement  at  the  centre  of  gravity  of  the  pile  cap  versus  the 

excitation  frequency  (first  mode). 


Table  8.2.  Computed  and  measured  resonant  frequencies  and  damping  ratios. 


Mode  of 
Excitation 

Computed 
Resonant  Frequency 
(Hz) 

Computed 
Damping  Ratio 

Calculated 
Damping  Ratio 

Vertical 

44.0 

46.5 

0.04 

0.05 

Translational 

6.67 

6.50 

0.06 

0.04 
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CHAPTER) 

CONCLUSIONS 

A  simplified  3-D  modd  of  the  response  of  a  soil  continuum  to  horizontal  earfiiquake 
shaking  has  been  devdoped  which  can  simulate  the  important  aspects  of  the  seismic  reqx>nse  with 
very  good  accuracy. 

The  model  has  been  formulated  in  terms  of  finite  dements  and  adiq;>ted  for  the  dyiumic 
analyris  of  piles  by  the  incorporation  of  beam  dements.  Solutions  can  be  obtained  for  both  dastic 
and  nonlinear  soil  response.  Nonlinear  response  is  modeUed  by  maintaining  compatibility  between 
shear  strains  and  effective  moduli  and  damping  throughout  the  dyiuunic  aiudysis. 

The  modified  3-D  modd  and  its  extension  to  dynamic  analysis  of  piles  was  the  origiiud 
conception  of  Guoxi  Wu.  (1993). 

The  model  has  been  validated  for  elastic  response  uring  existing  exact  elastic  solutions. 
The  soil  continuum  modd  alone  has  been  validated  using  Wood’s  (1973)  exact  solution  for 
dynamic  pressures  against  rigid  walls.  The  pile-soil  model  has  been  validated  by  comparing  pile 
impedances  for  single  piles  computed  by  the  model  ^h  those  computed  by  Kaynia  and  Kausd 
(1982)  using  iliU  3-D  continuum  equations.  Agreement  between  modd  solutions  and  the  more 
exact  solution  is  very  good. 

In  the  nonlinear  mode,  the  model  has  been  validated  for  angle  pile  response  uang  data 
fi'om  strong  shaking  tests  on  single  pile  foundations  conducted  on  the  centrifiige  at  the  California 
Institute  of  Technology  (Gohl,  1991).  The  important  aspects  of  accderation  and  moment  time 
histories  were  simulated  weU  by  the  model  and  the  distribution  of  peak  moments  along  the  pile 
were  within  6%  of  the  measured  moments. 
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The  model  simulated  successfully  the  response  of  a  full  scale  Franld  type  pile  to  forced 
vibration.  The  test  was  conducted  at  the  Univeraty  of  British  Columbia  Pile  Research  Fadlity  by 
Sy  and  Siu  (1992). 

The  computational  time  for  conducting  analyses  has  been  greatly  reduced.  Thus,  the  main 
objective  of  the  Phase  n  studies  has  been  achieved. 

It  now  remains  to  extend  the  model  to  pile  groups  and  to  dynamic  effective  stress  analysis. 
The  latter  feature  will  allow  consideration  of  the  effects  of  high  porewater  pressures  on  response. 
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